Quadrupolar gravitational radiation as a test-bed for /(/^-gravity 



Mariafelicia De Laurentis 1,2 , Salvatore Capozziello 1,2 
1 Dipartimento di Scienze Fisiche, Universita di Napoli "Federico II", 2 INFN Sezione di Napoli, 
Compl. Univ. di Monte S. Angelo, Edificio G, Via Cinthia, 1-80126, Napoli, Italy. 

(Dated: January 19, 2013) 

The debate concerning the viability of /(i?)-gravity as a natural extension of General Relativity 
could be realistically addressed by using results coming from binary pulsars like PSR 1913+16. 
To this end, we develop a quadrupolar approach to the gravitational radiation for a class of ana- 
lytic /(i?)-models. We show that experimental results are compatible with a consistent range of 
/(i?)-models. This means that /(i?)-gravity is not ruled out by the observations and gravitational 
radiation (in strong field regime) could be a test-bed for such theories. 
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I. INTRODUCTION 

The discovery of binary pulsars PSR 1913+16 by Hulse 
and Taylor in 1974 [l| opened a new testing ground for 
General Relativity (GR). In fact its continuous observa- 
tion by Taylor and coworkers 0, led to an impres- 
sively accurate tracking of the orbital motion of the bi- 
nary system. Before this discovery, the only available 
testing ground for GR was the Solar System where the 
gravitational field is slowly varying and represents only a 
very small deformation of a flat space-time. As a conse- 
quence, Solar System tests can only prove the weak-field 
limit of GR. By contrast, binary systems containing com- 
pact objects as neutron stars (NS) or black holes (BH) 
involve space-time domains where the gravitational field 
is strong. Indeed, the gravitational field on the surface 
(i.e. ~ 2 GM/c 2 R) of a NS is of order 0.4, which is close 
to the one of a BH ~ 2 GM/c R = 1) and much larger 
than the gravitational field on surfaces of Solar System 
bodies: ~ (2GM/c 2 R) Q ~ 10" 6 , (2GM/c 2 R) (B ~ 10" 9 . 
In addition, the high stability of pulsar clocks has made 
it possible to monitor the dynamics of its orbital motion 
down to a precision allowing one to measure the small 
(~ (v/c) 5 ) orbital effects linked to the propagation of 
the gravitational field at the velocity of light between the 
pulsar and its companion. The recent discoveries of the 
double binary pulsars 0, Q has renewed the interest in 
the use of binary pulsars as extremely relevant test-beds 
of gravity theories. This means that it is worth recon- 
sidering in detail, i.e. at its foundation, the problem of 
motion also in relation to the problem of generation and 
detection of gravitational waves (GWs). In other words, 
the motion of sources could give further signatures to 
GWs and then it has to be carefully reconsidered. 

The achieved sensitivity levels and theoretical develop- 
ments are leading toward a general picture of GW phe- 
nomena that was not possible in the previous pioneer- 
ing era. Experimentally, several GW ground-based laser 
interferometer detectors have been built in the United 
States (LIGO) Q, Europe (VIRGO and GEO) @, @| 
and Japan (TAMA) and are now taking data at de- 
signed sensitivities. A laser-interferometer space antenna 
(LISA) [l^l might fly within the next decade. As results, 



we can hope that the next decade will witness the direct 
detection of gravitational waves opening the fields of GW 
astronomy and cosmology. Theoretical studies have been 
developed in parallel to the experimental activity. In par- 
ticular, mechanisms for the production of GWs, both in 
astrophysics and in cosmology. Templates on binary in- 
spiral (see e.g. [TlT - fl4l | ) and robust search algorithms 
have been developed for GWs sources [l5[. Furthermore 
conceptual and technical problems, related to the pro- 
duction of GWs by self-gravitating systems (such as coa- 
lescing binaries) have not been fully solved. This status 
of art suggests to reconsider the problems of motion and 
generation of GWs also with respect to alternative the- 
ories of gravity which seem realistic approaches to face 
several problems in astrophysics and cosmology. In par- 
ticular /(ii)-gravity seem a viable semi-classical scheme 
to overcome shortcomings related to infrared and ultravi- 
olet behaviors of the gravitational field [l(|. These the- 
ories are based on corrections and enlargements of the 
Einstein GR. Besides fundamental physics motivations, 
they have acquired interest in cosmology due to the fact 
that they "naturally" exhibit inflationary behaviors able 
to overcome the shortcomings of Standard Cosmological 
Model (based on GR). The related cosmological mod- 
els seem realistic and, several times, capable of matching 
with the observations [T34H, From a genuine astro- 
physical viewpoint, these Extended Theories of Gravity 
(ETGs) [2(| do not urgently require to find out candi- 
dates for dark energy and dark matter at fundamental 
level (till now they have not been detected!). The ap- 
proach is very conservative taking into account only the 
"actually observed" ingredients (i.e. gravity, radiation 
and baryonic matter); it is in full agreement with the 
early spirit of GR which could not act in the same way 
at all scales (see [2(| for a comprehensive review). In fact, 
GR has been successfully probed in the weak-field limit 
(e.g. Solar System experiments) and also in this case 
there is room for alternative theories of gravity which 
are not at all ruled out, as discussed in several recent 
studies (22l - [24l |. In particular, it is possible to show that 
several /(i?)-models could satisfy both cosmological and 
Solar System tests [25|, [2(| , could be constrained as the 
scalar-tensor theories and could give rise to new effects 
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capable of explaining anomalies also at local scales (see 
for example [27| and references therein). 

In this paper we study the quadrupolar gravitational 
radiation in /(i?)-gravity using the "linearized theory". 
It consists in expanding the field equations around the 
flat Minkowski metric. The field equations then reduce 
to linear wave equations from which radiation can be cal- 
culated. GR predicts radiation that, at the lowest order, 
is proportional to the third derivative of the quadrupole 
momentum of the mass-energy distribution. It is a con- 
sequence of conservation equations that the first deriva- 
tive of the monopole momentum and the second deriva- 
tive of the dipole momentum are zero. This means that 
the gravitational radiation is first seen at the quadrupole 
term. The dipole effects depends on the difference of 
the self-gravitational binding energy per unit mass for 
two bodies and it is thus dependent also on the inter- 
nal structures of the objects. When the objects are in 
circular orbits, the time variation of the scalar field at 
each object, due to the motion of the other, is zero and 
the dipole contributions consequently drop out. Under 
these circumstances the dominant surviving terms are of 
quadrupole order. In /(i?)-gravity the situation is differ- 
ent due to the presence of further degrees of freedom of 
the gravitational field [29M3l| | . However, GR has to be 
fully recovered as soon as f(R) — > R. This "compatibil- 
ity" with GR could be a test-bed for these ETGs. Here 
we develop expressions for quadrupole gravitational radi- 
ation in /(i?)-gravity using the weak field technique and 
apply these results, to binary systems as, for example, the 
well known PSR 1913+16. In this way, it is straightfor- 
ward to compare the GR-predictions with those of ETGs. 
The outline of the paper is the following. In Sec. [Hi we 
briefly introduce the weak field limit and field equations 
of /(i?)-gravity. Sees. IIIII and ITVl are devoted to the cal- 
culation of the conservation laws. Finally the application 
to PSR 1913+16 is developed in Sec. [V] Conclusions are 
drawn in Sec. IVII 



II. FIELD EQUATIONS AND 
POST-MINCOWSKIAN LIMIT OF /(^-GRAVITY 

The post-Minkowskian limit of any theory of grav- 
ity arises when the regime of small field is considered 
without any prescription on the propagation of the field. 
This case has to be clearly distinguished with respect to 
the Newtonian limit which, differently, requires both the 
small velocity and the weak field approximations. Of- 
ten, in literature, such a distinction is not clearly re- 
marked and several cases of pathological analysis can be 
accounted. The post-Minkowskian limit of GR gives rise 
to massless gravitational waves. An analogous study can 
be pursued considering, instead of the Hilbert-Einstein 
Lagrangian linear in the Ricci scalar R, a general func- 
tion f(R) The only assumption that we are going to 
do is that f(R) is an analytic function. The gravitational 
action is then 
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f(R) + XC Tl 



(1) 
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where X = — - — is the coupling, £ rn is the standard 

c 4 

matter Lagrangian and g is the determinant of the met- 
ric 1 . The field equations, in metric formalism, read 2 



(2) 



3D/'(i?) + f'(R)R - 2f{R) = —T , (3) 



with T, 



the energy momentum ten- 



sor of matter (T is the trace), f'(R) = an d = 

;(T ;cr . We adopt a (+,—,—,—) signature, while the con- 
ventions for Ricci's tensor is R^ = R° \, a v and R a = 
Tp^ + ... for the Riemann tensor, where 



T af3 = [9aa.fi + ff/Ja.a ~ 9aP,a) , (4) 

are the Christoffel symbols of the g^ metric. Actually, in 
order to perform a post-Minkowskian limit of field equa- 
tions, one has to perturb Eqs. © on the Minkowski 
background rj^. In such a case the invariant metric ele- 
ment becomes 



ds 2 = g aT dx° dx T = (j] aT + h aT )dx cr dx T , (5) 

with hpu small {0(h) 2 <C 1). We assume that the f(R)- 
Lagrangian is analytic (i.e. Taylor expandable) in term 
of the Ricci scalar, which means that 3 



f(R) - E ^^(R Ro) n ^ fo + f' Q R + btlR 2 + - -(6) 



The flat-Minkowski background is recovered for R = 
i? ~0. 

Field equations @ , at the first order of approximation 
in term of the perturbation [3~ij , become : 



1 Here we indicates with "," partial derivative and with " ;" co- 
variant derivative with regard to all Greek indices run from 
0, 3 and Latin indices run from 1, 3; g is the determinant. 

2 All considerations are developed here in metric formalism. From 
now on we assume physical units G = c = 1. 

3 for convenience we will use for the following calculations, / in- 
stead of f(R) 
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where /„ = — 



R=0 



the gravitational theories that present also massive ones. 
A similar situation emerges also in the Newtonian limit: 
the Newton potential is recovered only as the weak field 
limit of GR. In general, Yukawa- like corrections, and then 
characteristic interaction lengths, are present (33|. The 



and □ = a' a that effective mass is 



R=0 



(S^)- 1 = and then f Q ' has 



is now the standard d'Alembert operator of flat space- 
time. From the zero-order of Eqs.([2]), one gets /(0) = 0, 
while T pv is fixed at zero-order in Eq.([7]) since, in this 
perturbation scheme, the first order on Minkowski space 
has to be connected with the zero order of the standard 
matter energy momentum tensor 4 . The explicit expres- 
sions of the Ricci tensor and scalar, at the first order in 
the metric perturbation, read 



d(!) _ ua _ lni _ lu 
n, pv — ll'(jj i ^ cr 2 V 2 .A 4 " 

Ri 1 ) = h' aT - Uh 



(8) 



with h = b? a . Eqs. (O can be written in a more suitable 

fa 

form by introducing the constant £ = —-jp, that is 

fo 

H(dl v^){Kr aT Uh) = ^Tf) . (9) 

By choosing the transformation hau = huu ~ ^Vuv an d 
the gauge condition „ = 0, one obtains that field 
equations and the trace equation, respectively, read 5 



(10) 



Uh + 3£,U 2 h 



.* T (o) 

fo 



It is worth noticing that solving the previous system of 
equations, we find wavelike solutions with massless and 
massive contributions (29l . l3ll . |32| . The presence of the 
massive term is a feature emerging from the higher-order 
terms in /(i?)-gravity. Specifically, it is related to the fact 
that £ 0, which is null in GR where f{R) = R. This 
means that massless states are a particular case among 



to be negative in order to have physically defined states. 
It is easy to see that massive modes are directly related 
to the non-trivial structure of the trace equation Eq.Q. 
In GR, the Ricci scalar is univocally fixed being R = 
in vacuum and R oc p in presence of matter, where p is 
the matter-energy density [5^, [3TL IsH ■ The task is now 
to evaluate the related energy-momentum tensors. 



III. ENERGY-MOMENTUM TENSORS 



Let us assume that the source T M „ is localized in a 
finite region. Outside this region T pv = 0. 

Then, as a consequence of Eqs. ([S} and gauge condition, 
we have 



Uh, 



0. 



outside the region. There are several ways to define the 
energy-momentum tensor of the gravitational field. One 
is to consider R^ on the left-hand side of Eq.([2]) con- 
sisting of a series of correction terms in R^l /. In the 

development of Eq.([7]), R^} is on the left-hand side. The 
remaining higher order terms, which so far have been ig- 
nored, could be brought to the right-hand side. If the 
source region gives rise to a flux of energy in the form 
of GWs, it must be represented by these higher order 
terms. This is the geometric approach (35|. The other 
approach is to use the standard field theoretical meth- 
ods. The geometric and the field theoretical approaches 
are complementary. Some aspects of GWs physics can be 
better understood from the former approach, some from 
the latter, and to study GWs from both vantage points 
results in a deeper overall understanding. We use the 
latter approach to calculate the stress-energy tensor of 
the gravitational field. So one can extending the formal- 
ism to more general theories and obtain this quantity by 
varying the gravitational Lagrangian. In GR, this quan- 
tity is a pseudo-tensor and is typically referred to as the 
Landau-Lifshitz energy-momentum tensor (36| . 
In the case of /(i?)-gravity, we have 



4 This formalism descends from the theoretical setting of Newto- 
nian mechanics which requires the appropriate scheme of approx- 
imation when obtained from a more general relativistic theory. 
This scheme coincides with a gravity theory analyzed at the first 
order of perturbation in the curved spacetime metric. 

5 The gauge transformation is ft' = ft M „ — f Mi „ — f„ iM when we 
perform a coordinate transformation as x'^ = x M + £ with 
0(C 2 )< 1. To obtain the gauge and the validity of the field 
equations for both perturbation h pv and h pv , the have to 
satisfy the harmonic condition = 0. 



5 d xy^gf(R) =5 Id xC(g^, g^ tP , g^ tPr7 ) 



4 / dC dC 2 dC 
d x\ dx- h o^- 



9g P a dg pa .\ 4 dg paM 



d^xJ—gHP°bg pa = 



The Euler-Lagrange equations are then 



(12) 
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dC 

dg P a 



dC 

dg P a,x 



a 2 dL 

K dg P aM 



0, 



(13) 



which coincide with the field Eqs. ([2]) in vacuum. Actu- 
ally, even in the case of more general theories, it is possi- 
ble to define an energy-momentum tensor that turns out 
to be defined as follows : 
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dC 



dC 



dg P a,xe. 



dg P cr,x dgpoM 



g P o,a 



(14) 



This quantity, together with the energy-momentum 
tensor of matter l) t „, satisfies a conservation law as re- 
quired by the Bianchi identities. In fact, in presence of 

matter, one has H pv = —T pvi and then 



= -^J—gTP°g p ^ a 



and, as a consequence, 



(15) 



It is important to stress that the definition of the energy- 
momentum tensor in GR and in /(i?)-gravity are differ- 
ent. This discrepancy is due to the presence, in the sec- 
ond case, of higher than second order differential terms 
that cannot be discarded by means of a boundary inte- 
gration as it is done in GR. We have noticed that the 
effective Lagrangian of GR turns out to be the symmet- 
ric part of the Ricci scalar since the second order terms, 
present in the definition of R , can be removed by means 
of integration by parts. 

On the other hand, an analytic /(i?)-Lagrangian can 
be recast, at linear order, as / ~ f^R + J-(R), where the 
function T satisfies the condition: lim^^o J 7 — > R 2 ■ As 
a consequence, one can rewrite the explicit expression of 

as: 



fX _ f'y-A 
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g P a,a 



OR 



d 9, 



p<?M 



~g P a^a 



J-"R.f— ■ — g P a,a — SuJ~- 



dg 



pcr,\(; 



(20) 



The general expression of the Ricci scalar, obtained by 
splitting its linear (R*) and quadratic (i?) parts once a 
perturbed metric ([5]) is considered, is 



-g(t*+XT*)l x = 0, 



(16) 



that is the conservation law given by the Bianchi iden- 
tities. We can now write the expression of the energy- 
momentum tensor in term of the gravity action f(R) 
and its derivatives: 



t x a = r 



OR 



dg P a,x 
OR 



OR 



dgp<?M 



g P (7,a 



dgptr,X£ 



dR 

g P cr,ia \ - f"R,£ q 9p<r,a ~ S a f , 



cgpuM 



(17) 



It is worth noticing that t„ is a non-covariant quantity in 
GR while its generalization, in fourth order gravity, turns 
out to satisfy the covariance prescription of standard ten- 
sors (see also On the other hand, such an expression 
reduces to the Landau-Lifshitz energy-momentum tensor 
of GR as soon as f(R) = R, that is 



l ... 



1 fdC Gl 



^xdg, 



pcr,\ 



-g P o 



5*£ 



GB 



(18) 



where the GR Lagrangian has been considered in its ef- 
fective form, i.e. the symmetric part of the Ricci tensor, 
which effectively leads to the equations of motion, that 
is 



\ tier pu nv an) 



(19) 



R = a ftv (T p -T p ) + a pu (T p T a - T a T p ) = 
= R*+R, (21) 

(notice that £gr = — V — gR)- I n the case of GR ia| GR > 
the Landau-Lifshitz tensor presents a first non-vanishing 
term at order h? . A similar result can be obtained in the 
case of /(-R)-gravity. In fact, taking into account Eq. (|20p . 
one obtains that, at the lower order, reads : 



r 



X _ fi.X 
a\h 2 — J0 L a\ GR 



-fo'R* 



dR* 



dg P a,xi 



g P a,a 



dR* 
dg P a, 



po\A£ 



-9pcr,£c 



og P aM 2 



fo^a 



| OR + f<> 

dR 



R 



-d e [R 



dg P a,xt 



( dR* 1 ; 



gp<?,i 



(22) 



Considering the perturbed metric ([3]), we have R* ~ 
R^\ where i?W is defined as in (|5J|. In terms of h and 
r], we get 



OR" 



rjP^^S — rj^rj 1 " 7 



(23) 
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Clearly, the first significant term in Eq. ([22)) is of second 
order in the perturbation expansion. We can now write 
the expression of the energy-momentum tensor explicitly 
in term of the perturbation h; it is 



K ~ /o^| GR +/o'{(^V- D/l ) 



+ ^ x a {h p \ p(J -Uh) 



h p(y h xs - 



+h p \ pa x h. a + h x ^ a Dh^ - Dh' x h a } . 



(24) 



Considering the tilded perturbation metric hp V , the more 
compact form 



\h\Uh - \h a Uh x - h x a , a uh°- 



(Dh) 2 6 x 
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is achieved. 

As matter of facts, the energy-momentum tensor of 
the gravitational field, which expresses the energy trans- 
port during the propagation, has a natural generalization 
in the case of /(i?)-gravity. We have adopted here the 
Landau-Lifshitz definition but other approaches can be 
taken into account (38j. The general definition of t x , ob- 
tained above, consists of a sum of a GR contribution plus 
a term coming from /(ii)-gravity : 



fo^a 



fX 



(25) 



However, as soon as f(R) = R, we obtains t x = t x , GR . 
As a final remark, it is worth noticing that massive modes 
of gravitational field come out from t X y since Oh can be 
considered an effective scalar field moving in a potential: 
t x , in this case, represents a transport tensor. 

The expression for the gravitational tensor t x can be 
simplified by doing approximations valid far from the 
source region. Far from the source hp V will be, functions 
of a single scalar variable t' 



where 



t' = t-r, 



r 2 = XiX 1 . 



(26) 



(27) 



Such a scalar can be constructed from the vector x p by 
forming 



with 



t' = k x x x , 



— k — 1 , k[ — X\ , 



(28) 
(29) 



r 



(30) 



In the far field k\ can be considered as a constant vector, 
over a small region. That is, h pv will be almost plane. 

Any — variation or change of the unit vector x over points 
r 

in the region can be made arbitrarily small by choosing 
a region far from the source j35| . 

The functional dependency of solutions will be on the 
t' = t — r. This fact can be done by expressing all partials 
of h^v as 



hfiu,a — Qx a dt' — ^a^ct^Vi/ — k a hp U , (31) 



where 



hfiv — h pl/ (k\x ) — h pl/ (t ) . 



(32) 



here the dot indicate the derivative with respect to the 
dx x 

time and — — = cL. Since Tuu = outside the source 
ox a 

region, 



(33) 



in the far field [23, HJ. If Eq.fl3T]i is used in the first of 
these, we find 



Oh^ v — h^ U p^ p — | kph^ u J , p — kpk^hpy , (34) 
implying that 



k p k p = . 



(35) 



Therefore, from Eq. (|24[) , the energy-momentum tensor 
associated with the tensor part of the gravitational field 
is 

= fo (k X k a h p °h pa ) + f» {kpk^k^Kh^- 
~\~ kp k(j \iP k k(y\\ ■ ■ ~ k p kjj \\P S 'q, k p k(j \\P ~\~ 

+^kpk a h pa 5 x ah - k e k a h Xi nh + k x k a hah+ 

+ ^8 x kpk a h p,J nh - ^S x (Dh) 2 - kpk a k s h pa k a h x Z- 



-kpk a h pc T k x k a h + k a h x £Dh^ - Ulv x k a h 



(36) 



Now remember that 

h = 7]£\h x £ , h = 7 1(x h x Z , (37) 

and 

k x v^ = h > (38) 



(i 



and then 

k x k a h = k x k a ri fx h x t = k^k a h^ , (39) 
we can further simplify t x in the following way 

t« = fo (k x k a hP a h p ^ + ft (k p k a h"ak^k a h^+ 

■ k p kjj \\P k kfy h ~^2j^^ ^ c% ^ p a ■ ■ 
+k p k a k i h pa k a h^ + k p k lJ h p ' T k x k a h s , (40) 

we notice that the sixth and fifth terms of above equation 
are equal because 

k p k a k x k a h pa h = k p k a k x h pa V^ = 

= k p k a k i h p °h p ^ (41) 

the third and fourth are the same, and then t x reduces 
to 



ta ~ fo k a h p h ptJ ^ — /q ^k p k a h p 5 a k p k a h p 



' fo k a h p h p(y 
~fX (kpKh^^^kpKh^ 



(42) 



finally the energy momentum tensor assume the following 
form 

t x = f k x k a (hrhpa) - \f^ x (k p k a h p °) 2 . (43) 



GR 



f(R) 



To be more precise, the first term, depending on the 
choice of the constant /q, is the standard GR term, the 
second is the f(R) contribution. It is worth noticing that 
the order of derivative is increased of two degrees consis- 
tently to the fact that /(-R)-gravity is of fourth-order in 
the metric approach. 

dE 

Now we could compute the instantaneous — — using 

at 

the Eq.(|43|l as a basis. The effect on a binary system 
is is more evident if we consider the average flux of en- 
ergy away from the system. Suppose that the h pv waves 
can be represented by a discrete spectral representation. 
The periodicity T will be proportional to the inverse of 
the difference of the pair of frequency components in the 

dE 

wave. Therefore, we must to evaluate the average of —j^- 

over an interval equal to or greater than T j3f| [36[ . The 
instantaneous flux of energy through a surface of area 
r 2 dfl in the direction x 



dE 
~~dl 

and the average flux is 



r'dVLxH 



i,0i 



(44) 



(45) 



and then Eq. |4"31 becomes 

(ti) = (fak x k a (h p °h pa ) - l -f{8 x (k p k a h pc 

(46) 

Finally, we re- write (t x \ in terms of a function J pv 
defined to be 

J pu (x,t) ~ 4 / d 3^ ^(^-|^-^D f (47) 



(48) 



\x — X 
noting that 

^/iiv(^',£) t ^ii/(*^,^), 
and consequently 



fo P h pa J P J per h P hper J P 'J per (^^) 



to give 



(ta) — ( f()k X k a j P<J j pa nfo^ai^p^a) J P<T Jp(J 



(50) 



IV. MOMENTA AND CONSERVATION LAWS 

Let us now analyze the radiation in terms of multi- 
poles, that means to expand J pu in a Taylor series about 
tf = t-r. That is, 



J^(x,t) = - 
r 



+x / dTx'x' 



d*x'T py (x',t')+ 

+ 



Of 



1 f ^s2d 2 T p "{x,t)' 



where we have used 



and 



\x — x\ 



\x — x ~ r — x ■ x 



(51) 



(52) 



(53) 



for r >> \x'\. Let us define the following moments of the 
mass-energy distribution: 



(54) 



M(t) ~ / d 6 xT w {x,t), 



D k {t) ~ / d 3 xx k T 00 (x,t), 



Q ij (t) ~ / cfeVT 00 (ic,i)- 



(55) 



(56) 
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The conservation law becomes, in the weak field limit, 
T^, v =0, (57) 
and implies the relations [HI, HU 
1 d 2 



d 3 xP k (x, t) 



2dt 2 



d 3 xx j x k T 00 {x,t) = 



(58) 



J d 3 xT 0k (x,t) = §- t J d 3 xx k T ao (x,t) = 

= D k (t), (59) 



Of 



cfWT» u (£,t) = I d 3 xT ]k {x,t) 

\Q jk {t), 



(60) 



We use Eq. ([5T]) to write J 00 , J 0i and J y in terms of the 
momenta Eqs. (f54"|) - (|56p . First, from Eq.(|5ip we have 



J 00 (x 1 t)=A- 
r 



d 



d 3 x , T m {x l ,t')+ 



+x~ / d d x'x n T w (x',t')+ 
at 1 



1„ . d 2 



jXiXj-^ I d 3 x'x' l x^T 00 (x',t') + 



(61) 



For J uu , it is easy to obtain 



J uu (£,i) 



4 1 



M{t')+x i D i (t') + -x i x j Q ij (t') 



(62) 



For J° l we need only two terms of the expansion Eq. (f5Tj) 
in order to include terms up to the second momentum 



r\x,t) = 4 



d 3 x'T°\x',t')+ 







+x k -^ I d 3 x'x' k T 0l (x',t') 



(63) 



Eg.(|59"]) and Eq.® then give 



J°\x,t) =4 



D%t') + -x k Q ik (t') 



For J iJ only one term of Eq. (|5ip is required, being 

J ij (x,t) = 2-Q ij (t'). 
r 

The conservation law also implies that 
M = , D k = . 



(64) 



(65) 



(66) 



Furthermore, from Ea.(|6"2"l). Eq.(J53J|, and Eq.(|6"5]l we 

have 



J — 2 X i X n j 

r 



j oi = 2-x k Q lk , 
r 



J 13 = 2-Q 
r 



and consequently 



J 2 XiXj 

r 



J m = 2-x k > 
r 



JV = 2-Q ij 
r 



(67) 



(68) 



(69) 



(70) 



(71) 



(72) 



In order to evaluate Eq.(|50p. we require that 

j pa J pa = j°° JOO + 2J° l j 0l + j VJ Jij . (73) 

and 

J p ° J pa = J 00 Jqo + 2 J m Joi + J" Jij . (74) 
Pluggins Eas.([67 ]l -([72 11 into Eq.([75)l and {73]), we get 



J P<T J, 



2(x k Q %k ) (., :()■') ■ {6' r Q,) 



(75) 



In completely analogous way, we find 



J p ° 'J P a — r2 



X 2.X j 



-2 (■< ! . Q' ! ) ( : >-jQ") + (Q''Q _ 

"(76) 

When Eq.flTS]) and Eq.flTf)]) are put into Eq.([50|). we 
find 



XiXjQ 1 ^ - 2 (x k Q tk ^j (xjQ tj 



(t a ) — ( fok k a 



Q t3 Qi 



{^X-iXj Q 



2ir L Q' L ) (xj'Q*) + ('Q^Qij 



(77) 



8 



Using the result in Eq. (|4"5j) and integrating over all 
directions in order to compute the total average flux of 
energy due to the tensor wave, 



^\ = r 2 f dilx l (t M ) . 

"•t I (total) J 



(78) 



Note that 



i a (t m ) = x l k°k%..] = x l (-l)(-x%..] = {...] , (79) 



which simplify the evaluation of Eq. (|T8[) . Integration over 
all direction is accomplished readily with the help of 



and 



„ i „ ■ 47T 

dflx l x J = —Sij , 



4tt 



(80) 



<Mx l x : >x l x m = — {5 l0 5 hn - S u 5 jm ) . (81) 



The result is: 

dE 
~dt 

(total) 



GR 



/(H) 



(82) 



Precisely, for /p and /q -> §, Eq.(|8"2"]) becomes 



dE 



G 



(GR) 



(83) 



which is which is the well-known result of GR [28J, 136J . 
See also (4(| for the recovering of the correct GR-limit. 

An important remark is necessary at this point. 
Eq. (18211 can be written as 



dE 
~dt 

(total) 



Gfo 

60 



(84) 



where the massive mode contribution is evident. This 
means that this further term affects both the total energy 
release as well as the waveform. This could represent a 
further signature to investigate such theories in the GW 
strong-field regime. 



V. APPLICATION TO THE BINARY SYSTEMS: 
THE PSR 1913+16 CASE 

Observations coming from PSR 1913+16 can be used 
to fix bounds on f(R) parameters. This could be consider 



a new test to retain or exclude such theories beside the 
classical Solar System experiments adopted for GR (28| . 
For a binary system, we have to assume that the motion is 
Keplerian in the first approximation and we can average 
over orbital periods. Given a point mass m, Q lJ (t) is 



''(*) 



d 3 x.TVT 00 (x^) 



dx 1 dxj dx? m x?x? x 
xS (x 1 - x 1 ^)) 5 (x 2 - x 2 (t)) 8 (x 3 - x 3 {t)) , 

(85) 

where x is the integration variable and x is the position 
of the mass [28[. We define m as the pulsar mass, M 

GM 3 

the mass of the companion star, and fi = the 

(M + my 

reduced mass. This last definition will be used to account 
for the fact that m can be small with respect to M. Since 
the orbit is Keplerian, we can choose x 3 — being a 
planar motion. Then Eq. (|85p reduces to 



Q i] (0=0, for i and/or j = 3 , 



Q n (i) = mix 1 ®) 2 , Q 22 (t) = m(x 2 {t)f 



(86) 



(87) 



)12 



(t) = Q 12 (t)=m[x 1 (t))(x 2 (t)] , 



where the position in the orbital plane is a function of 
time (28[. We are goin g to work in a parametric repre- 
sentation of the motion |36l. [39|-[4l[ and then let us recast 
the variables as 



a (1 — e cos£ ) , t 



-(£ — e sm£) , 



x {£) = a (cos£ — e) , x (£) = a (l — e J 2 sin£ , 

(89) 

where r is the orbital radius, a, the semi-major axis of 
the orbit, e, the eccentricity, £ . the eccentricity anomaly. 
Over the whole orbit, £ ranges from to 2ir. We use 
£, rather than t, to locate the body in its orbit, and 
therefore we have to integrate over d£ 



For a Keplerian orbit T has the value 




(90) 



(91) 



Therefore, if f(t) = g(£(t)), we may write Eq. (|90f as 

(/(*)} = £- P 9(£)(l-ecos£)d£. (92) 
27T Jo 



9 



Q can be expressed in terms of the eccentric anomaly 
and then Eq. (|92|) can be used to compute the time aver- 
age over an orbital period 6 . We find that time derivative 
can be recast as 



d d£ d 2n d 
Jt = Hd£ = T {1 - £COs£r df 



From Eqs.(|86 



Q ij Q = (Q 11 



we can write 



2 (Q 1 



Q 



(93) 



(94) 



Let us consider the various orders of derivation. First, 
from Eq.JHZI) and Ea.(f89l we have 



Q n (£) = ma 2 (cos£~e) 2 . (95) 
Using Eq. (f9"3")) to compute the derivatives, we find 
o /27r\ sin£ (cos£ — e) 



n (f) 



-2ma 



T 



1 



: COs£ 



(96) 



in 



f y ) 



X 

x (2 cos 2 £ - e cos £ - e cos 3 £ + e 2 - l) 



and 



finally 



= —2ma 



(1 — e cos£) 5 sin£ x 



x (e cos 2 £ + 2e 2 cos £ - 4 cos £e 3 + 4e) 



(97) 



-2ma 2 ( ^ ) (! ~ ecos£) 7 x 



[(16c 



Be 2 + 4) cos 2£ + (8e 2 - 3e) cos £+ 



+3e(cos3£-4e(2e + l))] . 

Likewise, from Eq. ([57|) and Eq. (|89")) . we have 

Q 22 £ = ma 2 (l - e 2 ) sin 2 £ , 
which leads to the derivatives 



i'2it\ , 9 , sin£cos£ 
Q 2 \£) = 2ma 2 [-) (l - e 2 ) ^—^ , 



(98) 
(99) 



Note that we can rise/lower space indices without regard for sign 
changes because rf J = 5 lJ . 



'■(£) = 2ma 2 ( % 



l-e 2 



T J (l-ecos£) 3 ' 
(cos 2 £ — sin 2 £ — e cos 3 £ ) 



-22 



(£) = 2m,a 



2n 



TJ (l-ecos£) 
x sin £ (3e - 4 cos £ + e 2 £ 



f x 



(100) 



2 ' 2?r 



'(£) = 2ma 2 — 



1-e 2 



(1-ecosf) 7 
x [(22e 2 - 16)cos2£ + 41ecos£+ 
+ (e 2 (cos4£-39)-9ecos3£)] (101) 

Finally, from Eq.flSS]) and Eq.® we have 

Q 12 = Q 21 = ma 2 (l - e)3 sin f (cos f - e) , 



(102) 



whose derivatives are 
'2tt 



Q 



(1 - e)*(l - ecosf)" 1 x 
x(2cos 2 £ — ecosf — f) , 



ni2 



Q 



x (2e cos 2 £ - 4 cos £ + 2e) , 



<r = 



I J (1 — e)5(l — ecosf) 5 x 

x (e 2 cos 2 £ + 3ecos£ + ecos 3 £— 
-3e 2 -4cos 2 £ + 2) , 



mfl2 ( ) (1 - (1 - ecos£)~ 1 



(103) 



T 



x sin 



£ [(I5e 



4) cos£ + (3e 3 + 6e) cos2£- 



+27e 3 + e 2 cos3£ + 18e] 



(104) 



When results from Eas.([97]l. ([TU0]| . and (flQg]) . together 
with Eq.(|9Tf for T, are used in Eq.([94"]t. one finds 



4m 



2 ^ [8(l-e) + 6 2 sin 2 g] 
(1 — ecos£) 6 



(105) 



10 



and 



Q V Q\ 



2m 



2 M 



1 



a 8 (1 - ecosf) 



14 



x2(e 2 - l) 2 [41ecosf - 9 - 39e 2 + 

+ (22e 2 - 16) cos 28 + e cos AS] 2 + 

+(1 - e) [I8e - 27e 3 + (4 + 15e 2 ) cos£+ 

+(6e + 3e 3 ) cos 28 + e 2 cos 38] 2 sin 2 £ + 

+2 [(8e 2 - 3e) cos£ + (4 + 8e 2 + 16e 3 ) cos£+ 

2ecos3f - 16e 3 -8e 2 ] 2 

(106) 

Substituting Eq. ([T05]) into Eq.j92]) and averaging, we 
have 



The above results apply for the motion of a body of 
mass m in a Keplerian orbit about a second body of mass 
M. Therefore, we can evaluate the overall loss rate due 
to the motion of both bodies. Let the subscript 1 denote 
the position of the pulsar m and 2 that of the companion 
M and, as above, let the coordinate origin be at the 
barycenter (421 . [43[ . This condition gives 



and then 



mx\ + Mxl = . 



M 



(III) 



(112) 



The averall momentum Q lJ for the system consisting of 
both m and M is 



4 ™2M 3 8 (1 - e) + e 2 sin 2 f 
(1 — e cos£ ) 5 



dS. 



Q' 



mx\x{ + Mx\x> 2 = — (m + M)x\x{ , (113) 



(107) 



also for Eq. (|106p we obtain 



where we have used Eq. (|112[ ) to express the momentum 
in terms of the motion of m. Averaging for the binary 
system, we obtain 



Q tr Q\ 



m 2 ^ 



\13 



lo (1 — ecosf) 
x2(e 2 - l) 2 [41ecos£ - 9 - 39e 2 + 

+ (22e 2 - 16)cos2£ + ecos4£] 2 + 
+ (1 - e) [I8e - 27e 3 + (4 + 15e 2 ) cos£+ 
+ (6e + 3e 3 ) cos 28 + e 2 cos 3£] 2 sin 2 8 + 
+2 [(8e 2 - 3e) cos 8 + (4 + 8e 2 + 16e 3 ) cosf - 



Q"Q, 



G 6 mrM 1 25e 4 + 196e 2 + 64 
2a 5 (m + M) 4 (i - e 2)5 



(114) 



d£ 



We do not measure — directly but the change of orbital 
at 
d8 

period T induced by — . To this end, we remember that 
the semi-major axis of the orbit is (36l . l39l . |421 - [44| 



2ecos3£"- 16e 3 - 8e 2 l 2 



, m + M 



(115) 



The first term of Eq. (|107[) is evaluates using 



^ \there it has to be recalled that a is the semi-major axis 
of the pulsar orbit. The total energy E of a Keplerian 
binary system is then 



d£ 



o (l-ecos£) 5 (l- e 2) 

where P 4 (x) = - (35a; 4 - 30x 2 + 3) . 
8 

dS tt 3e 4 + 24e 2 



E 



GmM 
2a' 



GmM 2 
2a(m + M) 



from which 



o (l-ecos£) 5 8 (1-e 2 )^ 



The complete evaluation of Eq. (|107[) is 



GmM 2 1 
'2(m + M)E ' 



(116) 



(117) 



(109) 



The orbital period T can be related to the energy E by 
combining Eq. (|9Tj) and Eq. (|117| . The result is 



Q° J Q ij ) = ^m 2 



1 2 /i 3 25e 4 + 196e 2 - 
-m —r — 

2 a 5 (1-e 2 )^ 



64 



T 



-ttGE? 



3 M 3 



2(m + M) 



(118) 



(110) 



and for the Eq. (|108[) we do not have an analytical solu- 
tion of the integral but, only a numerical result that will 
be insert in the following equations. 



By taking the time derivative and Eq. to restore the 
parameter a, we find that 



dT 
It 



(m + Mf 



a 5 dE 
G 3 M 7 It 



(119) 
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Let us now use the published numerical values for the 
specific example of PSR 1913+16 to numerically evalu- 
ate the above equations . The results will be included 

,'dE\ „ . . dT 

into Eq. (|82[) to evaluate 



dt 



from which — can be 
dt 



estimated using Eq. (|119l) . We use the values from Taylor 
et al. P,[3| for PSR 1913+16 reported in Table I. 



PSR 1913+16 


chacteristic features 


pulsar mass 


m = 1.39M 


companion mass 


M = 1.44M Q 


inclination angle 


sin i = 0.81 


orbit semimajor axis 


a = 8.67 x 10 10 cm 


eccentricity 


e = 0.617155 


gravitational constant 


G = 6.67 x 10~ 8 dyn cm 2 g~ 2 


speed of light 


c= 2.99 x lO^cmsec" 1 



Table I: Values from Taylor et al. for PSR 1913+16 
First we find, from Eq.fpII]), that for PSR 1913+16 



T = 2.21 x 1(T 



Q 



2.78 x 10 y 



dE 



g cm 2 J dt 



92 gem* 



see-' 



urx =1.29xl0 96 ^ Cm 



(120) 



fl211 



(122) 



Then, the averaged total radiation rate for the binary 
system in f(R)- gravity is found from Eg. ([52")) , that is 7 



34 f/n gem 



, = [1.32 x 10 31 /o - 6-10 x 10 34 /o] 

Using this value in Eq. (|120p . we find 
f f{R) = 2.92 x 10- 13 ./o - 1.34 x IQ- 9 .^' (^) , 



(123) 



(124) 



as we can see from the above equation the orbital period 

depends strongly on the choice theory. Now, given the 
4 

value of /q, (i.e. /q = -), we determine the value of /q 
3 

that falls within the limits observed by Hulse and Taylor 
P, 0| ■ We remember that they predicted an upper and 
lower limit in the observation of the orbital period that 



7 Note that for dimensional reasons we have to restore the factor 
of c 5 to the denominator of . 

dt 



is about 3.8 x 10" 12 < t < 2.6 x 10~ 12 and the limit for 
GR is 



T G r ~ 3.36 x 10" iZ — 



sec\ 



sec/ 



(125) 



We immediately recover GR limit putting /q = — and 

3 

ftf = 0. In Fig. 1 is shown a plot of f from (fT24|) 
for PSR 1913+16 as a function of /q parameter. In 




-dT(f(R)) 
GR limit 

Lower limit set by HT 
- Upper limit set by HT 



-1.4 



-1.2 -1 

x to' 3 



Figure 1: Orbital decay rate for PSR 1913+16 in f(R)- 
gravity. Upper limit set by Taylor et al. in dashed line. GR 
limit 3.36 x 10~ 12 in dotted line and the lower limit set by 
Taylor et al. in dashdot line. While in solid line is plotted 



T 



/(«)■ 



Fig. 1, the observational limits on T are indicated to- 
gether with the GR limiting value (|125j) . The range 
-2.63 x 10" 3 < fH < -2.25 x 10~ 3 well fits with these 
observational limits [l|, Q. In other words, we can con- 
clude that /(i?)-gravity is not excluded by the Hulse and 
Taylor observations on binary pulsar. On the other hand, 
such observations contribute to fix the range of viability 
of such theories. 



VI. CONCLUDING REMARKS 

In this paper, we developed the post-Minkowskian 
limit of analytic /(i?)-gravity models in the Jordan frame 
to calculate the gravitational radiation emitted by a bi- 
nary system. One of the results is that the quadrupole- 
radiation, in /(i?)-gravity and in GR, occurs indepen- 
dently of the detailed internal structure of the stellar 
bodies. It depends on the masses of the two bodies, on 
the orbital parameters and on the details of the gravita- 
tional theory. Further massive modes emerge and they 
are directly related to the analytic parameters of f(R)- 
gravity, that is the coefficients /g and /g of the Taylor 
expansion. This fact is relevant since it does not depend 
on specific /(i?)-models but it is a general feature. 

As a consequence, the theoretical quadrupole radia- 
tion rate, calculated according to the theory, can be con- 
fronted to binary system observations to fix the param- 
eters of the theory. Specifically, the radiation rate is a 
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function of f' and /g'. As we can see from Fig. 1 or, 
equivalently from Eq. (|124[) . the predicted range of the 
time derivative of the orbital period for PSR 1913+16 
is compatible with the observational uncertainty estab- 
lished by Hulse and Taylor. [45[. This means that obser- 
vations can fix the parameters of the theory. These re- 
sults pose interesting problems related to the strict valid- 
ity of GR. It seems that it works very well at local scales 
(Solar System) where effects of further gravitational de- 
grees of freedom cannot be detected. As soon as one is 
investigating larger scales, as those of galaxies, clusters 
of galaxies, etc., further corrections can be introduced in 
order to explain both astrophysical large-scale dynamics 
(4(| and cosmic evolution [471 . |48| . Alternatively, huge 
amounts of dark matter and dark energy have to be in- 
voked to explain the phenomenology, but, up today there 
are no final evidences for these new constituents at fun- 
damental level. What we have shown is that the Hulse 



and Taylor experiment, beside confirming GR, does not 
exclude Extended Theories of Gravity [20| including GR 
MS i\ particular case. 

Furthermore, the fact that, up to now, only massless 
gravitational waves have been investigated could be a 
shortcoming preventing the possibility to find out other 
forms of gravitational waves. Tests in this sense could 
come out, for example, from the stochastic background 
of gravitational waves where massive modes could play a 
crucial role in the cosmic background spectrume [49J, l50l] . 
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